AO-A036  71«  TEXAS  TECH  UNIV  LUBBOCK  OEPT  OF  ELECTRICAL  EN6INEERINB  F/6  12/1 


A NEM  CHARACTERIZATION  OF  THE  NYOUIST  STABILITY  CRITERION* (U) 

AU6  76  R OECARLO*  R SAEKS  AF-AF0SR-2631-7A 

UNCLASSIFIED  AFOSR-TR-77-0069  NL 


MICROCOPY  RtSOLUHON  UST  CHART 

NA' ■ NAi  BIIPIAII  uf  !'«-•  A 


MILWAUKEE  C? 


AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH  (AFSC) 
NOTICE  OF  TRANSMITTAL  TO  DEC 

Xhi3  tochnical  report  bar.  bejii  reviewed  and  i3 
approved  for  public  release  IA<V  Ai'R  190-lH  (7b). 
Distribution  is  liulimited. 

A.  D.  BLOSE 

Technical  Information  Officer 


A NEW  CHARACTERIZATION  OF  THE  NYQUIST  STABILITY  CRITERION’ 


R.  DECARLO  & R.  SAEKS 


ABSTRACT 


The  usual  proof  of  the  Nyquist  Theorem  depends  heavily  on  the  agrument  principle. 

The  argument  supplies  unneeded  information  in  that  it  counts  the  number  of  encircle- 
ments of  "-1".  Stability  of  a system  requires  an  encirclement  or  a no-encircle- 
ment test.  Using  homotopy  theory,  this  paper  offers  a more  intuitive  approach.  We 
believe  this  approach  v/ill  lead  to  practical  generalizations.  For  example,  systems 
characterized  by  several  complex  variables  such  as  multi  deminsional  digital  filters. 


I.  INTRODUCTION 

This  paper  introduces  a characterization  of  the 
Nyquist  criterion  using  homotopy  theory,  a branch 
of  algebraic  topology.  The  authors  emphasize  the 
intuition  and  motivation  for  this  approach.  The 
hope  is  to  aid  interested  readers  to  further  extend 
and  apply  these  ideas.  In  this  vein,  proofs  are 
omitted  so  as  to  simplify  the  presentation.  Details 
can  be  found  in  the  references.  With  this  phil- 
osophy in  mind,  let  us  define  the  type  of  system 
we  will  be  discussing. 

As  illustrated  in  Figure  1,  let  g(s)  be  a rational 
function  in  the  complex  variable  s (bounded  at  s = 
„♦*)  representing  the  open  loop  gain  of  scalar 
single  loop  feedback  system.  The  closed  loop 
system  has  transfer  function  h(s)  = g(s)/[l+g(s)]. 
The  closed  loop  system  is  stable  if  and  only  if 
all  poles  oj  h(s)  are  in  the  open  left  half  plane 
denoted  by  !Z_  (where  t will  denote  the  entire  com- 
plex pl.ane). 

The  Nyquist  Criterion  states  that  the  closed  loop 
system  is  stable  if  and  only  if  the  Nyquist  plot 
of  g{s)  (i.c.  the  image  of  the  Nyquist  contour 
under  the  map  g(-))  does  not  encircle  nor  pass 
through  ”-l".  If  the  Nyquist  plot  passes  through 
"-1"  there  is  a pole  on  the  imaginary  axis;  if  the 


**This  boundedness  condition  can  bo  dispensed 
with  & is  added  only  to  ease  the  exposition. 


Nyquist  plot  encircles  "-1",  there  is  a pole  in  the 
gpen  right  half  plane,  which  we  will  denote  by 
2^  will  denote  the  closed  right  half  plane).  ■ 
The  following  section  constructs  the  required  mach- 
inery of  homotopy  theory 


II. 


Figure  1 

MATHEMATICAL  PRELIMINARIES  & BACK- 
GROUND 

Basic  to  homotopy  theory  is  the  concept  of  a path. 

A path  or  a curve  in  the  complex  plane  is  a contin- 
uous function  of  bounded  variation  (2)  y :[o,l]«?. 

Y is  a closed  path  if  y(0)  ' y(1)-  y is  a simple 
closed  path  if  y is  a closed  path  and  has  no  self 
intersections.  The  image  of  I = [0,1]  under  y is 
called  the  trace  of  y and  is  denoted  by  Iy). 
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Figure  2 

Two  closed  curves  and  homotopic  in  £ if 

there  exists  a continuous  function  rilxl^?  such 
that: 

(a)  r(s,0)  = y^(s)  0 < s<  1 

(b)  r(s,l)  = y^(s)  0 < s<  1 

(c)  r(C,t)  = r(l,t)  0 < t<  1 

Intuitively,  y^  is  homotopic  to  y^  if  one  can  con- 
tinuously deform  y^  into  y-j.  Moreover,  it  is  easily 
shown  that  the  homotopy  relation  is  an  equivalence 
relation.  (4)  (5) 

Another  important  property  of  a closed  curve  is 
its  index  or  degree.  The  index  (2)  is  of  closed 
curve,  with  respect  to  a point  "a"  not  in  {yl  is; 


dz 


t 


i 


This  integral  measures  the  net  increase  in  angle 
that  the  ray  r of  Figure  3 accumulates  as  its  tip 
traverses  the  trace  of  y. 


Intuition  for  the  approach  stems  in  part  from  the 
observation  that  n(y;-l)  = 0 if  and  only  if  y is 
homotopic  to  a point  in  £-  {-11  (cf.  prop.  5.4, 
ref.  2).  We  will  henceforth  refer  to  such  a y as 
being  homotopical ly  trivial.  Conversely,  y en- 
circles "-I"  if  and  only  if  y cannot  be  continuous- 
ly deformed  to  a point  in  £-  {-!).  These  ideas 
appear  to  indicate  that  the  Nyquist  encirclement 
condition  is  fundamentally  a homotopy  concept. 

The  intuition  is  further  reenforced  when  one  form- 
ulates the  Nyquist  criterion  on  the  Riemann  sur- 
face (2)  (8)  associated  with  a map,  f(s).  Assum- 

ing f(s)  is  analytic  on  and  bounded  at  s=-»,  the 
image  of  simply  connected  regions  in  are  simply 
connected  in  £.  To  illustrate  the  point,  let 


Figure  4-a  be  the  image  of  the  right  half  plane 
under  f(s).  The  region  is  not  simply  connected. 
Figure  4-b  shows  the  “same  region"  as  it  might 
appear  on  an  appropriate  Riemann  surface.  Here  the 
region  is  simply  connected. 


The  boundary  of  the  regions  depicted  in  Figure  4 
are  the  Nyquist  plots  of  f(s)  in  I and  on  the  Rie- 
mann surface.  On  the  Riemann  surface  the  Nyquist 
test  becomes  an  obvious  triviality.  In  £ it  is 
mathematically  more  delicate. 

Our  setting  uses  homotopy  theory,  a branch  of  alge- 
braic topology,  to  establish  a topologically  in- 
variant relationship  between  a metric  space,  X,  and 
an  algebraic  group  called  the  fundamental  group  of 
X,  denoted  by  n(X).  The  relationship  is  topologic- 
ally invariant  in  that  homeomorphic  spaces  have 
isomorphic  fundamental  groups. 

Specifi'-ally,  the  fundamental  group  is  a set  of 
equivalence  classes  of  closed  curves.  Each  equiva- 
lence class  consists  of  a set  of  curves  homotopi- 
cally  equivalent.  The  group  operation  is  "concate- 
nation" of  curves. 

Fo  example,  the  fundamental  group  of  £ consists  of 
one  element,  i^,  the  identity,  since  all  closed 
curves  are  homotopic  to  zero.  If  X = £ - {-1), 
then  ii(X)  has  a countable  number  of  elements:  ij^ 
(the  identity)  equal  to  the  equivalence  class  of 
all  closed  curves  not  encircling  "-1"  and  the  re- 
maining elements,  (n  = 1 , 2,  3 . . . ) consisting 
of  the  equivalence  class  of  all  closed  curves  en- 
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circling  "-1",  n times.  Moreover,  concatenated 
with  is  equal  to  the  element 
Now  let  X and  Y be  metric  spaces.  Let  f;  X->-Y  be 
locally  homeomorphic.  In  particular,  assume  that 
for  each  point  y in  Y there  exists  an  open  neigh- 
borhood G of  y such  that  each  connected  component 
of  f'^(G)  is  homeomorphic  to  G under  the  map  f. 
Under  this  condition  X is  said  to  be  a covering 
space  of  Y.  (2)  (4)  Also  let  ti(X)  and  t!{Y)  be  the 
fundamental  groups  associated  with  X and  Y respect- 
ively. With  these  assumptions,  f effects  a group 
isomorphism  (i.e.  a one  to  one  into  mapping  pre- 
serving group  operations)  between  n(X)  and  a 
subgroup  of  i7(y)  as  in  the  following  diagram  (4) 

(5) 

X F 'TT(X') 

1 


a critical  value  of  g,  but  g(a)=  0 with  g’(a)  i 0. 
Finally,  define  T = ftjt.  = f\-l),  i=l,...,nl. 
Note  also  that  since  f is  a rational  function,  P, 

Q and  T are  finite  sets.  Define  X = {!  - {PuQuT) 
and  define  Y = f (X) . 

lerma  1 : Under  the  above  hypothesis,  X is  a cover- 

ing space  of  Y.  This  leads  to  the  following  corol- 
lary. 

Corollary:  The  fundamental  group  n(X)  of  X is  iso- 

rorphic  to  a subgroup  N of  n(Y). 

This  corollary  says  that  a closed  curve  in  X is 
homotopically  trivial. 

III.  THE  SCALAR  CASE 

Let  g(s)  be  as  described  in  the  introduction.  Appro- 
priately define  the  sets  P,  Q,  and  T and  the  spaces 
X and  Y so  that  X is  a covering  space  of  Y.  Also 
as  per  reference  (10)  and  Figure  6,  construct  the 
ugly  Nyquist  contour,  and  the  usual  Nyquist 
Contour,  r,  where  r;l  ? o {<»). 


% 

McTT(y) 


Figure  5 

F is  the  functor  which  establishes  the  relation- 
ship between  a topological  space  and  its  fundamen- 
tal group.  Finally  let  us  distinguish  between  a 
critical  point  and  a critical  value.  A point 
in  ? is  a critical  point  of  a differentiable  func- 
tion f if  f'(Zg)  =0.  A critical  value  of  f is 
any  point  w=  f(Zj,)  whenever  z^  is  a critical 
Point. 

Now  suppose  f:?  is  a rational  function  whose 

set  of  polos  is  P = IP^ Pj^l.  Let  Q=(qj 

q^j^l  be  the  set  of  all  points  in  0 such  that  f(q.) 
is  a critical  value  of  f.  Note  that  there  may  be 
q.'s  which  are  not  critical  points.  To  see  this 

' O 

consider  g(z)  =z  (z-a).  g (0)  =0  implies  "0"  is 


lu;  /*' 


X indicates  a point  of  P;  S indicates  a point  of 

Q 


Figure  6-(b) 

X indicate?  a point  of  P;  0 indicates  a point  of  Q 
Lctunia  2:  Under  the  above  assu:iiptions  on  g and 

fi(s)  is  stable  if  and  only  if  the  path  goXj^ 
does  not  encircle  "-1".  (10) 

At  this  point  we  must  establish  this  lemma's  con- 
nection with  the  classical  Myquist  criterion.  To 
this  end  we  compare  the  information  of  the  Nyquist 
plot,  gor  with  the  "ugly"  Nyquist  plot, 

Lemma  3:  Let  n be  the  number  of  poles  of  g in 

then 


r 

"if 

1 (z-l)-'dz  = 

2TT 

gop 

fii'se  three  lemmas  give  rise  to  the  following 
theorem. 

Theorem  1:  Let  g(s)  be  as  above.  Then  h(s)  is 

stabel  if  and  only  if  the  Nyquist  plot  of  g(s) 
does  not  pass  through  "-1"  and  encircles  "-1" 
exactly  ^n  times  where  n is  the  number  of  poles  of 
g(s)  in 


IV.  MATRIX  CASE 

Let  the  entries  of  an  nxn  matrix  G(s)  be  rational 
functions  in  the  complex  variable  s.  Suppose  G(s) 
characterizes  the  open  loop  gain  of  the  single  loop 
feedback  system  of  Figure  7. 


Figure  7 

x(s)  and  y(s)  are  n vectors  whose  entries  are  also 
rational  functions  of  s which  represent  the  input 
and  output  of  the  system  respectively. 

This  article  assumes  each  entry  of  G(s)  is  bounded 
at  X = «>.;  Thus  G(s)  as  a mapping,  G(-):? 
is  analytic  on  % except  at  a finite  number  of 
points,  the  poles  of  its  entries. 

For  Figure  7 to  be  well  defined  we  require  that  det 
[I+G(s)]=0  Thus  there  exists  a closed  loop  convolu- 
tion operator,  H,  such  that  y = U*x.  Moreover  the 
Laplace  transform  of  H,  H(s)  satisfied 

H(s)  = G(s)[I+G(s)r'' 

For  this  system  to  be  stable,  H(s)  must  have  all 
its  poles  in  and  have  all  its  entries  bounded 
at  s = «■. 

Under  the  assumptions  on  G(s),  the  following  factor- 
ization is  valid: 

G(s)  = N(s)D‘\s) 

where  M(s)  and  0(s)  are  right  co-prime,  polynomial 
matrices  in  s with  det[D(s)]70.  Moreover  s^  is 
a pole  of  G(s)  if  and  only  if  it  is  a zero  of  det 
[D(s)].  (9) 

Desoer  and  Schulman  (3)  have  shown  that  the  close 
loop  operator  H is  stable  if  and  only  if  det[N(s)+ 
D(s)]^0  for  s in  and  det[I+G('«)]70.  Using 
this  fact,  we  state  and  prove  the  following: 

Theorem  2:  H is  stable  if  and  only  if  (1)  the 
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Nyquist  plot  of  det[N(s)+D(s)]  does  not  encircle 
nor  pass  through  "O'',  and  (2)  det[I+G(“)]7‘0.  (10) 

Observe  that  if  one  assumes  the  open  Igop  gain  to 
be  stable  (i.e.  G(s)  has  all  poles  in  i^)  then 
det[I+G(s)]  in  the  above  theorem.  This  follows 
since  for  all  s in  (5^,  det[N(s)+D(s)]  = det[I+G(s)] 
det[D(s)]  with  det[D(s)]^0.  Thus  in  det[N(s)+ 
D(s)]  has  a zero  if  and  only  if  det[l+G(s)]  has  a 
zero. 

Finally,  it  is  worthwhile  to  point  out  the  rela- 
tionship between  the  above  formulated  multivariable 
Nyquist  criterion  and  that  formulated  by  Barman 
and  Katznelson.  For  this  purpose  we  let  Aj(iw); 
j=l,  ....  n;  denote  the  n eigenvalues  of  G(iw).  In 
general  parameterization  of  these  function  by  iw  is 
not  uniquely  determined  but  one  can  always  form- 
ulate such  d function.  Moreover  these  functions 
are  piecewise  analytic  and  can  be  concatonated  to- 
gether in  such  a way  as  to  form  a closed  curve 
which  BaiTan  and  Katznelson  term  the  Nyquist  plot 
of  G(s). 

Now,  Since 

n 

det[l  + G(iw)]  = n [1  + (iw)] 

and  tlie  degrees  of  a product  is  the  sum  of  the 
degrees  of  the  individual  factors  and  also  equals 
the  degree  of  the  cancetonation  of  the  factors, 
the  degree  of  the  Barman  and  Katznelson  plot  with 
respect  to  "-1"  coincides  with  the  degree  of  our 
plot  with  respect  to  "0".  As  such,  even  though 
the  two  plo'.s  are  different  their  degrees  coincide 
and  (lence  either  can  be  used  for  a stability  test. 
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